function P=LM_man_step(G,P,k)
P_new=P;
while 1
    P_new=s_random(G,P,k);
    if inpoly(P_new.pos,G.boundary.nodes),
        break;
    end;
end;
P=P_new;
